Phase Behavior of Bent-Core Molecules 



(N 
O 

o 

(N 
X> 

IX, 



o 

£ 

■4— » 

c3 



^3 

G 

O 

o 



(N 
O 
<N 
O 
-i— > 



i 

c 

o 
o 



X 



Yves Lansac, Prabal K. Maiti, Noel A. Clark, and Matthew A. Glaser 
Condensed Matter Laboratory, Department of Physics, and 
Ferroelectric Liquid Crystal Materials Research Center, 
University of Colorado, Boulder, CO 80309, USA 
(Dated:) 

Recently, a new class of smectic liquid crystal phases (SmCP phases) characterized by the sponta- 
neous formation of macroscopic chiral domains from achiral bent-core molecules has been discovered. 
We have carried out Monte Carlo simulations of a minimal hard spherocylinder dimer model to in- 
vestigate the role of excluded volume interations in determining the phase behavior of bent-core 
materials and to probe the molecular origins of polar and chiral symmetry breaking. We present the 
phase diagram as a function of pressure or density and dimer opening angle ip. With decreasing ip, 
a transition from a nonpolar to a polar smectic phase is observed near ip = 167°, and the nematic 
phase becomes thermodynamically unstable for ip < 135°. No chiral smectic or biaxial nematic 
phases were found. 

PACS numbers: 61.30.-v, 61.30.Cz, 64.70.Md 



Molecular chirality plays an important role in the sci- 
ence of liquid crystals (LCs), leading to cholesteric LCs 
blue phases [0, ferroelectric || and antiferroelectric 
[I| smectic phases, and twist grain boundary phases ||. 
In all of these examples, chirality is an intrinsic property 
built into the chemical structure of the LC molecules. Re- 
cently, a new class of smectic LC phases (SmCP phases) 
characterized by the spontaneous formation of macro- 
scopic chiral layers from achiral molecules has been dis- 
covered [|| |t|. The molecules comprising these phases 
have 'bow' or 'banana' shaped cores. The resulting 
phases exhibit two spontaneous symmetry-breaking in- 
stabilities: polar molecular orientational ordering within 
the layer plane, and molecular tilt, which together pro- 
duce chiral layers with a handedness that depends on the 
direction of the tilt relative to the polar axis. Very large 
second order nonlinear optical (NLO) coefficients have 
been measured in the ferroelectric state of such materi- 
als, bearing some promising applications in NLO devices 

From a theoretical point of view, the relationship of 
phase behavior to the specific bent-core molecular shape 
is of fundamental interest. In this paper, we investigate a 
minimal excluded volume model of bent-core mesogens, 
focusing on the molecular origin of polar and/or chiral 
symmetry breaking. Of particular interest is the coupling 
between polar and chiral symmetry breaking. In all bent- 
core materials studied to date, polar symmetry breaking 
is accompanied by chiral symmetry breaking induced by 
molecular tilt. This empirical fact raises the question 
whether there is a fundamental connection between po- 
larity and chirality in molecular fluids. Another empirical 
observation is that bent-core materials exhibiting SmCP 
phases generally do not exhibit nematic phases, although 
two exceptions have recently been reported |l0|, [fl| . One 
objective of this study is to establish the molecular shape 
requirements for the occurence of the nematic phase in 
bent-core materials. Finally, we explore the possibility 



of biaxial nematic ordering in bent-core materials, moti- 
vated by recent experimental indications |ll|, |l2| ]. 

Hard core models are particularly appealing due to 
their simplicity and relative ease of computation, both 
in simulation and theory. In particular, hard sphero- 
cylinders have been widely studied as simple models for 
conventional LCs |b| |l4j . This model exhibits rich phase 
behavior including isotropic, nematic, smectic, columnar 
and solid phases, the phase transitions being driven by 
the competition of two main cntropic contributions, the 
orientational entropy favouring the isotropic phase and 
the positional entropy favouring ordered phases, as shown 
in the forties by Onsager in the limit of infinitely thin 
rods U|. 

To capture the main characteristics of the collective 
behavior of bent-core molecules, we extend the sphero- 
cylinder model by introducing a hard-core dimer formed 
by two interdigitated hard-core spherocylinders sharing 
one spherical end cap (see inset in Figure |l|) . This is an 
ideal model system to consider due to the relatively small 
parameter space. There are three parameters: two geo- 
metrical parameters, namely the length-to-breadth ratio 
L/D, where L and D are respectively the length and 
the diameter of each spherocylinder, and the opening 
angle ip between the two spherocylinder axes; and one 
thermodynamic parameter, the reduced pressure P* de- 
fined as P* = PPv or, equivalently, the reduced den- 
sity p* defined as p* — pv D . Here, v Q is the volume of 
the equivalent straight hard spherocylinder (tp = 180°), 
v = ttD 3 /6 + nL'D 2 /A with L' = 2L. In all the simula- 
tions presented below, we consider a single value for the 
length-to-breadth ratio, namely L/D — 5. This ratio has 
been chosen in order to roughly represent the equivalent 
geometrical envelope described by a realistic bow shaped 
molecule, with fully extended aliphatic tails. 

Using a similar model, Camp et al. |l6| have found 
a nematic phase and a smectic A phase of bent-core 
molecules roughly half as long as the ones considered 
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here. However, no systematic study of the phase behavior 
was carried out. Recently, a polar smectic A phase and a 
chiral crystal phase have been observed for a 'polybead' 
model of bent-core molecules with an opening angle of 
ip = 140° @. 

To investigate the phase behavior of our model as a 
function of the pressure and of the opening angle ip, we 
perform NPT Monte Carlo (MC) simulations, with peri- 
odic boundary conditions, on a system of N = 400 bent- 
core molecules. For each opening angle, the system is 
initially prepared at high pressure in the crystal phase 
corresponding to the highest number density (antipolar 
crystal, see Figure ||). The equation of state P*{p*) is 
determined by measuring the density p* as the pressure 
is decreased incrementally from the crystalline state. At 
each state point (ip, P*), 2 x 10 5 MC cycles are used for 
equilibration and 1 x 10 6 MC cycles are used for produc- 
tion of the results. 

The location of the phase boundaries is determined 
from the equation of state P*(p*), and the nature of 
the coexisting phases is investigated through the compu- 
tation of various order parameters. In-layer crystalline 
order is probed with the translational order parame- 
ters pGk defined as pck = jt J2j=i exp(iG k ■ Tj) where 
Gi, G2 and G3 are the reciprocal basis vector of a de- 
formed hexagonal lattice, r\,- is the position of the cen- 
ter of mass of the molecule j, and M is the number of 
molecules in a given layer. The smectic order is mea- 
sured by the layer translational order parameter pu de- 



fined as 
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where d is the layer spacing and z is the layer nor- 
mal. Polar orientational order is detected using the 
polar order parameter m defined as m = 5 J2j=i ™j 
where m 3 is a unit vector contained in the plane of the 
molecule and passing through the apex of the molecule 
JL8| . Finally, molecular orientational order is deter- 
mined from the largest eigenvalue of the second-rank 
tensorial orientational order parameter Q a/ 3, defined as 
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\5 a p), where a,fi = x,y,z, 



and rij is a unit vector parallel to the molecular end-to- 
end vector for molecule j. 

The (ip,P*) and (ip,p*) phase diagrams are presented 
in Figure Rich phase behavior is found, with isotropic 
(I), nematic (N), smectic A (SmA), polar smectic A 
(SmAP), columnar (Col), crystal (X) and polar crystal 
(XP) phases. Configurations from the isotropic, nematic, 
polar smectic and polar crystalline phases are shown in 
Figure [|, for an opening angle of tp = 165°. 

The nematic phase is stable for opening angles larger 
than tp — 135°. With decreasing opening angle, the re- 
gion of stability of the nematic phase decreases, vanishing 
for opening angles smaller than ~ 135°, leading to an (I, 
N, SmAP) triple point near ip = 135°. It is interesting to 
note that for dimers half as long as ours, small opening 
angles seem to destabilize the smectic phase rather than 
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FIG. 1: Phase diagram of spherocylinder dimers (inset) 
with aspect ratio L/D — 5 as a function of opening angle ip 
and reduced pressure P* (top) and reduced density p* (bot- 
tom). All two-phase regions are shaded. The following phases 
are present: isotropic liquid (I), nematic (N), polar smectic 
A (SmAP), smectic A (SmA), columnar (Col), polar crystal 
(XP) and crystal (X). 



the nematic phase Jig ]. 

The vast majority of bent-core materials exhibiting liq- 
uid crystal behavior, and in particular SmCP phases, 
have an opening angle between 120° and 135°, and do 
not exhibit any nematic phase. Quite interestingly, two 
classes of bent-core compounds having an opening an- 
gle between 134° and 148° JjlJ [H| exhibit both smectic 
and nematic phases. These observations, in good qualita- 
tive agreement with the predictions of the model, tend to 
confirm the hypothesis that excluded volume interactions 
play a central role in the behavior of such materials. 

The existence of a biaxial nematic phase remains an 
elusive possibility in thermotropic LCs. Due to their 
bent-core geometry, banana molecules are good candi- 
dates to investigate such phenomena. Recent experi- 
ments suggest that the nematic phase exhibited by two 
classes of bent-core material might be biaxial Jill p"2[. 
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FIG. 2: Final configurations from Monte Carlo simulations 
of N = 400 bent-core molecules with opening angle ip = 165° 
as a function of pressure. From left to right: isotropic phase 
(P* = 1), nematic phase (P* — 5), polar smectic A (P* — 
10), and polar crystal (P* — 15). 

However, the nematic phase presented by our hard-core 
model does not exhibit any biaxiality. It is likely that 
the reported biaxiality is due to more subtle interactions 
and/or to the presence of the flexible tails. A biaxial 
nematic phase has been reported for hard spherocylin- 
der dimers twice as long as the ones used in the present 
study, but no transition from/to an isotropic liquid has 
been reported Q . Using a simplification of the Onsager 
second- virial treatment and bifurcation analysis, Teixeira 
et al. have found a biaxial nematic phase in the limit of 
very long bent-core molecules Ipof . 

Because straight spherocylinders do not exhibit any 
polar smectic ordering, it is expected that our model 
should exhibit a transition from nonpolar smectic (SmA) 
to polar smectic (SmAP). This transition occurs for an 
opening angle between 167° and 169°, and is associated 
with two triple points, a (SmA, SmAP, N) triple point 
near ip = 167° and a (SmA, SmAP, XP) triple point near 
ip = 169°. Very recently, the first group of bent-core 
molecules exhibiting both SmCP and a SmA (as well as 
SmC and nematic phases in one case) have been synthe- 
sized 0, The opening angles measured by NMR 
techniques for the region of appearence of the smectic A 
remains in the range of 132° to 145° but no polar or- 
der has been detected (^2), lowering the upper limit of 
stability of a polar smectic phase with respect to our 
predictions. 

Figure |^ displays the evolution of the order parame- 
ters for opening angles ip — 157.5° and ip — 172.5°. In 
the former polar smectic phase characterized by 

a high value of both the smectic and polar order param- 
eters is present while in the later case the appearance of 
a smectic phase from the higher density crystal phase is 
accompanied by a jump to zero of both the crystal order 
parameters and the polar order parameter. The smec- 




reduced pressure 



FIG. 3: Evolution of the squared-magnitude of order param- 
eters as a function of reduced pressure for an opening angles of 
ip = 157.5° (top) and ip = 172.5° (bottom) showing, respec- 
tively the phase sequences XP-SmAP-N-I and XP-SmA-N-I 
as a function of decreasing pressure. The following order pa- 
rameters are plotted: (A, V, <) solid-liquid order parameters 
Pi, p2, pz; (□) polar order parameter m; (O) smectic order 
parameter pu; (O) the largest eigenvalue of the nematic order 
parameter Q a p. 



tic phases (both polar and nonpolar) do not exhibit any 
tilt. Neverthless, our results tend to confirm the hypoth- 
esis that the polar ordering is related to the distinctive 
bent-core shape of the molecules but is not strongly re- 
lated to the molecular tilt. Thus, polarity does not imply 
chirality. 

A transition between a polar crystalline phase and a 
narrow nonpolar crystalline phase (i.e a rotator phase) 
is also present. This rotator phase is stable for opening 
angles larger than ip = 172.5°, and is characterized by 
a (SmA, X, XP) triple point around ip = 172.5°. Quite 
interestingly, the rotator phase competes with a colum- 
nar phase for opening angles larger than ip — 174° and 
smaller than ip — 180°. This narrow columnar phase is 
characterized by a significant crystal order parameters 
but a negligible magnitude of the smectic order param- 
eter. Since no evidence of such a phase was found for 
straight spherocylinders, a slightly bent molecular shape 
seems to favor the appearence and the stabilization of 
a columnar phase. Due to the rather unexpected ap- 
pearence of the columnar phase, we performed additional 
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studies of this region of the phase diagram using helical 
periodic boundary conditions p3| for N = 400 and a di- 
rect 'quench' from a crystalline state to the middle of the 
columnar phase for tp = 176° with helical periodic condi- 
tions for a larger system (N = 1600). A columnar phase 
was observed in both studies. However, we feel that free 
energy computations are needed in order to assess the 
relative thermodynamic stability of the columnar phase, 
the nonpolar smectic phase, and the rotator phase. 

Insights into the shape of the phase boundaries can 
be gain by supposing, to a first approximation, that the 
partition function of the system can be decomposed into 
a product of a positional and orientational contributions, 
in which case the entropy is the sum of an orientational 
entropy and a translational entropy. Competition be- 
tween different forms of entropy determines the stabil- 
ity of a given phase at a given density. In the limit 
of straight spherocylinders, the isotropic-nematic phase 
transition occurs when the gain in positional entropy S pos 
exceeds the loss of orientational entropy S orlcnt [fl5|| . A 
nematic-smectic phase transition occurs when the gain in 
translational entropy perpendicular to the long molecular 
axis S^° s exceeds the loss of positional entropy parallel 
to the long molecular axis S^, leading to the forma- 
tion of a stack of two-dimensional liquid layers. Simi- 
lar reasoning can be applied to bent-core molecules: in 
the range 134° < tp < 180°, the isotropic phase is more 
favourable at smaller opening angles. As the cores be- 
come more bent, the gain in positional entropy asso- 
ciated with nematic ordering is reduced. The nematic 
phase range is then reduced, eventually disappearing for 
tp < 134°. The shape of the nematic-SmAP boundary 
(i.e., for 134° < tp < 167°) can be qualitatively under- 
stood in the same way by noticing that the positional 
entropy parallel to the long molecular axis S^ os is larger 
for larger opening angles than for smaller ones, stabiliz- 
ing the nematic phase for larger opening angles. This 
trend is reversed for the nematic-SmA transition (i.e., 
for tp > 167°) because the absence of polar order leads to 
jamming, reducing the translational entropy perpendicu- 
lar to the long molecular axis S , *^ ans for decreasing open- 
ing angles. This effect is responsible for the enhanced 
stability of the nematic phase for decreasing opening an- 
gles. 

The spherocylinder dimer model exhibits rich phase 
behavior, including polar and nonpolar crystal, colum- 
nar, polar and nonpolar smectic, nematic and isotropic 
phases. In particular the existence and range of stability 
of the nematic phase is in good agreement with the be- 
havior of the new class of bent-core molecules, while the 
stability of the nonpolar smectic A phase is in qualita- 
tive agreement with experiments. In addition our model 
predicts the existence and stability of a polar smectic A 
phase for tp < 167°. Such phases remain to be found. No 
tilted phases are exhibited by the model, and a simple 



extension of the present model, taking into account the 
important steric role played by the flexible LC tails in 
the formation of tilted smectic phases, is currently un- 
der investigation. Our model indicates that there is no 
intrinsic coupling between polar symmetry breaking and 
chiral symmetry breaking, and that the later is not di- 
rectly related to the bent-core shape of the molecules. 

This work was supported by NSF MRSEC Grant DMR 
98-09555. 
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